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PRIMITIVE POINTS IN LATTICE POLYGONS 


IMRE BArANY, GREG MARTIN, ERIG NASLUND, AND SINAI ROBINS 


Abstract. Let A be a star-shaped polygon in the plane, with rational vertices, con¬ 
taining the origin. The number of primitive lattice points in the dilate tA is asymptoti¬ 
cally Area(<A) as t —>■ oo. We show that the error term is both (ti/log log t) and 
0(t(logt)^/^(loglogt)^'^^). Both bounds extend (to the above class of polygons) known re¬ 
sults for the isosceles right triangle, which appear in the literature as bounds for the error 
term in the summatory function for Euler’s 4’{n). 


1. Introduction 

One of the fundamental problems in discrete geometry is to estimate the number of lattice 
points contained in a polygon. In this paper we concern ourselves with the set of primitive 
lattice points 

(1.1) P = {{m,n) e gcd(m,n) = 1}, 

also known as lattice points visible from the origin. It is a classical result the number of 
primitive lattice points in a “reasonable” region in is approximately times the area of 
the region. We shall be interested in the family {tA} of dilates of a hxed polygon A, for 
which we dehne the error term 

(1.2) E_^{t) = #(tAnP)-;rArea(fA) = ^{tAn¥) --f^Area(A). 

The fact that ^(tAnF) ~ Area(A), or equivalently that E^^it) = o(f^), was likely used 
as far back as Minkowski (see [12, page 998] or [2, Theorem 459]). Stronger upper bounds 
for Ej\^{t) are known; we state the following result, which will be justihed in the next section, 
as a benchmark for comparison. 

Proposition 1.1. If A is a polygon, then 

Ej^{t) <^tlogt 

for t > 2, where the implicit constant may depend on A. 

The purpose of this paper is to improve this upper bound for any rational polygons (a polygon 
all of whose vertices have both coordinates rational), and show that it cannot be improved 
too much more by providing a strong result for the error term. 

It is instructive to consider the specihc example A = A, where A is the isosceles right 
triangle with vertices (0, 0), (1, 0), and (1,1) and thus area Then ^(fAflP) is the number 
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of primitive points in the dilate tA = {(x, ?/) G : 0 < ?/ < x < t}, that is, 

(1.3) #(tAnP)= 1 = 1+ E 

0<m<t 0<n<m l<m<t 

gcd(m,72)=l 

(where the extra 1 comes from the fact that gcd(0,1) = 1). It is well known that this 
summatory function of the Euler 0-function is asymptotic to so we define 

(1.4) EA(t) = ^ -h 1 = #(tA n P) - 

l<m<t 

Proposition 1.1 implies the estimate E^(t) tlogt, which (when phrased in terms of the 
summatory function of the Euler phi function) is a classical result of Mertens [5]. The best 
unconditional upper bound for this error term E^(t) is due to Walfisz [13, page 144, eq. (3)]. 
Using methods related to exponential sums, he showed that 

(1.5) E^it) < f(logf)+^(loglogt)+^. 

As a consequence of our work, we can extend this bound from the isosceles right triangle A 
to all rational polygons. 

Theorem 1.2. Let .4, C be a rational polygon. Then 

Ej,{t) <^t{\ogtf/\\og\ogtf/\ 
where the implicit constant may depend on A. 


Prior to Walfisz’s result, Chowla and Pillai [10] showed that these upper bounds cannot 
be improved too much by establishing the lower bound 


lim sup 

t^OO 


\EA{t)\ 
t log log log t 


> 0 , 


that is, E^(t) = (flog log log f). Later, Erdos and Shapiro [1] obtained a slightly weaker 
quantitative lower bound but showed that E/^{t) oscillates in sign infinitely often. Both 
results were improved by Montgomery [6]; 


Proposition 1.3. E^it) = f2±(f\/loglogf); in other words, 

^ ^ j V ■ e E^it) 


lim sup _ 

t^oo t^/\og\ogt 


> 0 and lim inf — ^ < 0. 

t^oo t^/\og\ogt 


We use the term origin-star-shaped (star-shaped with respect to the origin) to refer to any 
domain B for which XB C B for all A G [0,1]; note that in particular, any origin-star-shaped 
domain contains the origin. The main result of this paper generalizes this lower bound of 
Montgomery to this class of polygons: 


Theorem 1.4. Let M C be a rational origin-star-shaped polygon. Then 

E_A{t) = fi±(fVlog logf), 

where the implicit constant may depend on A. 
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In the process of applying an adaptation of Montgomery’s argument to the general error 
term Ey\^{t) for the primitive point counting function for the lattice polygon A, we found 
ourselves needing to establish the following “error term independence” result concerning the 
error term for the summatory function of 0(n), which was dehned in equation (1.4); this 
result may be of independent interest. 

Theorem 1.5. For any positive rational numbers Ci,... ,Ck and /i,..., fk, 

CiE^ifix) H-h CkE/^{fkx) = ff±(x0og log x), 

where the implied constant may depend upon the Cj and fj. 

In other words, there can be no “magic cancellation” among the terms E^{fjx) that makes 
the oscillation of the sum signihcantly smaller than that of an individual term. Furthermore 
we conjecture that oscillations of the same size exist even when the Cj are allowed to be 
negative (where we require the /j to be distinct to avoid trivial cancellations). 

Conjecture 1.6. Theorem 1.5 holds for any rational numbers Ci and any distinct positive 
rational numbers fi. 

This would imply a stronger version of Theorem 1.4 that holds for any rational polygon 
not necessarily star-shaped or containing the origin: 

Conjecture 1.7. Let M C be a rational polygon. Then 

Ejy{t) = Q±{t^yloglogt), 

where the implicit constant may depend on A. 

The rest of the paper is divided into two sections. In the next section, we show (Theo¬ 
rem 2.6) that the error term Ej[{t) may be rewritten as a linear combination of dilates of the 
totient error function E/^(t), thereby establishing Theorem 1.2 and reducing Theorem 1.4 to 
Theorem 1.5. We then establish this latter theorem in the hnal section. 

2. Decomposing the error term 

We begin by giving a proof of Proposition 1.1, not only for the sake of completeness, 
but also because the structure of the argument foreshadows how we will approach the main 
result of this section, namely Theorem 2.6 below. We use the term pointed triangle to mean 
a triangle which has the origin as a vertex. 

Proof of Proposition 1.1. We begin by quoting a reasonably precise estimate [3, special case 
of Theorem 2.1] for the number of primitive points inside a domain B\ if B is convex and 
contains the origin, then 

6 

(2.1) #(i3nP)- -ATea{B) max{l,a;loga;}, 

where u is the diameter of B. In particular, if C is a convex polygon containing the origin 
and B = tC is a dilate, then the diameter of S is a constant multiple of t and so Ec{t) 
max{l,flogt} (with the implicit constant depending on C). 

However, any polygon A can be partitioned into sums and differences of hnitely many 
pointed triangles, simply by triangulating A and noting that any triangle can be writ¬ 
ten as the signed sum of three pointed triangles. Let denote these triangles, and 
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Sj e {—1,1} the sign of Cj, that is whether it is added or substracted, so that Area(y4) = 
Yl]j=i Area(Cj). Even though these triangles have sides in common, the number of double- 
counted lattice points on the f-dilation of each side is at most 0{t), and therefore the above 
bound (valid since pointed triangles are certainly convex polygons containing the origin) 
implies 

6 / ^ \ 6 
EA{t) = hi P) - - Area(f^) = ( fl P) -|- 0{t) j- - Area(t^) 

\ j=i ' 

^ / g \ g 

= ^ e J — Area(tCj) + {t) j + 0{t) - - Area(M) 

j=i ^ 

n 

= ^ ejEc^ (f) + 0(f) < f log t 
i=i 

for f > 2, as desired. □ 

Remark. In [3], the bound (2.1) is stated without the hypothesis that B contain the origin; 
however, this hypothesis is actually necessary. One can easily construct, using the Chinese re¬ 
mainder theorem, a square B of diameter cj containing no primitive lattice points whatsoever. 
The area of such a square is a constant times cj^, and therefore \^{Bn¥) — ^ Area(i3) | S> cj^, 
which is incompatible with the claimed bound cjlogcj. 

The error in the proof of [3, Theorem 2.1] comes when applying Lemma 2.3, which requires 
CJ > 1, to a term of the form X](A/fc)i /(^^)) when k > u, the diameter of A/k is less than 
1, and so Lemma 2.3 cannot be applied. However, if A is a set containing the origin, then 
these sets A/k are sets with diameter less than 1 that contain the origin, hence contain no 
primitive lattice points (indeed, no lattice points at all other than the origin). Therefore the 
sums over k in the proof of [3, Theorem 2.1] can be truncated at fc < cj, and the rest of the 
argument goes through thereafter. 


The remainder of this section is dedicated to proving Theorem 2.6, which asserts that the 
error term E^it) for any rational polygon A may be rewritten in terms of the totient error 
function E/^{t). We begin with a detailed investigation of this latter function. 


Definition 2.1. Throughout this section, we employ the notation [xj for the greatest integer 
not exceeding x and {x} = x — [xj for the fractional part of x. We also employ the sawtooth 
function dehned as 

B{x) = {x} - ^ 

for all real numbers x; this function is equal to the hrst periodic Bernoulli polynomial .Bi(x) 
except at integer arguments. 


First we establish, for the sake of completeness, some elementary, well-known, and useful 
estimates for sums involving the Mobius p-function, which satishes the key identity 


( 2 . 2 ) 




1, if n = 1, 

0, otherwise. 


Lemma 2.2. For any real numbers 2 < x < y: 
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E 

d<x 


fi{d) 6 


la) > —Tj— — + O 

71 ^ 


(b) 

d<y 


(c) 


E 

d<x 


d 


= 1 ; 
< 1 . 


Proof. We begin by noting that 


E 

d>x 


fi{d) 


d^ 


^ 1 dt 

- = 

d>x 


1 1 

_iV X — 1 X 


since x >2. Thus part (a) follows from the classical evaluation 


p{d) 

^ d^ 

d=l 


6 


C(2) vr2 


Next, since = J2j<z have 


d<y 


X 


j<x d\j 2<j<x 


d<y j<^ld j<x d<y 

d\j 


by the identity (2.2), establishing part (b). Lastly, note that 


X 




d<x 


d<x 


+ 


— 1 + y^/r(d)| 

d<x ^ 


by part (b). The right hand side is bounded in absolute value by 1 + a: <C x, establishing 
part (c). □ 


We are now ready to give a more precise formula for the totient error term E^it). 
Proposition 2.3. For any real numbers T >t >2, 

N / 
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Proof. Starting from equation (1.3) and the key identity (2.2), we have 
#(tAnP) = i+ 1 


l<m<t l<n<m 
gcd(m,n)=l 

m<t n<m d|gcd(m,n) 
gcd(m,72)=l 

d<t m<t n<m d<t k<tld i<k 

d\m d\n 


d<T k<t/d d<T 

2 . 


+ 


1 




d<T 


+ 


by Lemma 2.2(b). In the remaining sum, we expand — (2 ~ {^})^ obtain 


d<T d<T '' d<T '' ■' 


1 

2 


tV6 

= 2b+o' = 


by Lemma 2.2(a) and a trivial bound. Therefore 
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EA{t) = #(tA n P) - —r Area(A) 


TT^ 


A+2 _ + pjd) J _t 1 

7r2 ^ d \ d j 

/^(^) 


+ 0(T)-— 


TT^ 


d<T 


d 


B{PP-]+om 




d<T 


by Lemma 2.2(c). 


□ 


Counting lattice points in the simplest way (each with weight 1) is problematic in our 
present context. First, we will be decomposing polygons into unions of triangles that share 
sides, and so double-counting lattice points on these shared boundaries would become an 
issue; second, any error term as large as the perimeter for counting lattice points in a triangle 
would result in an unacceptably large error term in the inclusion-exclusion method we use 
to detect primitive lattice points. Consquently, we employ a more convenient weighting in 
our lattice-point counting, namely the solid angle sum of a polygon. 
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Definition 2.4. For a point p G let B{p]r) = {q : d{p — q) < r} denote the ball of radius 
r around p. Let A denote Lebesgue measure on and define for any polygon A 


uja{p) = lini 

r^O 


\{B{p] r) n A) 
\{B{p-r)) 


It is easy to work out that 


ua{p) = < 


fo> 

1 , 

1 / 2 , 

l^p/27r, 


if p ^ Al, 

if p is in the interior of A, 

if p lies in the interior of an edge of A, 

if p is a vertex of A, 


where Op is the angle at the vertex p. 

Using this function, we dehne the solid angle sum of the f-dilate of the rational polygon 
A to be 


We also dehne the corresponding sum over primitive lattice points only: 

^A{t) = ^UJtAip)- 

peP 

The utility of using the solid angle weighting is that both of these functions are additive, in 
the sense that A^ub(^) = ^A{t) + ^b(^) and P^ub(^) = IPyl(^) + IPb(^) for any polygons A 
and B with disjoint interiors. 

Proposition 2.5. Let A he a rational polygon. There exist a real number C, a positive 
integer k, and rational numbers Ci,..., c^, /i,..., /^ where fi>0 (all depending on A) such 
that 

k 

(2.3) AA{t) = Area(Al)f^ + Ct — t''^^CjB{fjt) + 0(1), 

i=i 

where the implicit constant may depend upon A. Furthermore, when A is a origin-star-shaped 
Polygon the Cj may he taken to all be positive. 

Proof. Le Quang and Robins [4] gave a precise formula for the solid angle sum of any rational 
pointed triangle T: 

k(T) 

(2.4) Ar{t) = Area(r)f2 + C{T)t - t ^ cfiT)B{ffiT)t) + 0(1) 

i=i 

where 0(T), k{T), Cj(T), fj(T) are constants depending on T with k{T) a positive integer, 
and the Cj{T) and fj{T) positive rational numbers. But the rational polygon A can be 
partitioned as a signed sum of a hnite number of rational pointed triangles (as in the proof 
of Proposition 1.1). Simply summing the formula (2.4) over these hnitely many triangles 
yields the desired formula (2.3), where O is the sum of the 0(T), and { fA is the union of 
the {/j(r)}, and so on. □ 

The following theorem results from a careful examination of P^(f), which can be related 
to AaA) using the Mobius function as in the proof of Proposition 2.3. 
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Theorem 2.6. Let A he a rational polygon. There exist a positive integer k and rational 
numbers ri,... ,rk, fi,..., fk where /* > 0 such that 

k 

(2.5) E^(t) = 5^r,EA(/,t) + 0(t), 

i=i 

where the implicit constant may depend upon A. Furthermore, if A is a origin-star-shaped 
polygon, then the rj may all be taken to he positive. 

Remark. Theorem 1.2 follows immediately from Theorem 2.6 in light of the known upper 
bound (1.5). In addition, Theorem 1.4 follows immediately from the combination of Theo¬ 
rem 2.6 and Theorem 1.5; the latter theorem is proved in the next section. 

Proof. We commence by excluding the (non-primitive) origin and using the Mobius identity 
(2.2) to write 

P^(t) = fi{d). 

(m,,n)sZ2\(0,0) d\{m,n) 

Interchanging the order of summation (valid since in reality there are only finitely many 
nonzero terms) and rescaling (which preserves solid angles), 


OO 


P. 4 (() = 

E 

u}tA{m,n) 

d=l 

(m,n)EZ^\(0,0) 



d|m, d\n 


CXD 

= 

E ‘ 


d=l 

(3;,y)sZ2\(0,0) 


OO 

= 5^+(c^) 

d=l 

(■43 - 

-a;t_ 4 ( 0 , 0 ) 


by the definition of We may truncate the outer sum at any real number T >t diam(.4,) 
is allowed because for larger values of d, the diameter of ^A is less than 1, and hence j^A 
(which contains the origin) cannot contain any other lattice points. Proposition 2.5 now 
implies 

P.4(t) = 5^/i(c?)(Area(.A)^ + C'^-^^c,-5('/,^') +0(1)) 

d<T ^ “ j=l ^ 

d<T d<T 

j=l d<T ^ ' ^d<T ' 

= Area(.4)t2(4 + o(^ij) +0(T), 

3 —1 d^T 






Using Proposition 2.3, we may rewrite the above in terms of E/^(t) and obtain 

(2.6) ¥At) = Area(^)4t' + jE^.t) + 0(T) 

^ i=i E 

for any T > f max{l, diam(^), /i, ..., fk}. 

The number of integer points on the boundary of tA is 0{t) since tA has finitely many 
sides, each of which has length 0{t). Consequently, 

P.4(t) -#(MnP) 

and therefore equation (2.6) implies 
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EaW = #{tA n P) — Area(^)— 

71^ 

= P^(f) + 0(f) - Area(^)-f2 = + 0(T). 

^ j=i E 

Upon setting T = f max{l, diam(^), /i, ..., /fc} and rj = Cj/fj for each I < j < k, the 
theorem follows. □ 


3. Linear combinations of Ea 

In this section we prove Theorem 1.5, showing that positive rational linear combinations 
of scaled copies of the totient error function Ea have oscillations as large as those known 
for Ea itself. We begin by recalling some of the components of Montgomery’s argument [6] 
establishing the oscillations of Ea, after which we describe the strategy that led to our 
modifications. 


Definition 3.1. Define 


Roix) = 


(f){n) 6 


n<x 


n 71 


-—X. 


Montgomery [6, Theorem 1] showed that the totient error function is closely connected to 
the above weighted error: 

Lemma 3.2 (Montgomery). Ea{x) = xRq{x) + 0(a:exp(—c\/logx)). 

Definition 3.3. Define 

K{q,a) = 


d\q 


d 


dr 


where the sawtooth function B{x) was defined in Definition 2.1, and 

-1 


C{q,a) = K{q,a)-^Yl ( ^ “ 


p\g 




Lemma 3.4. If b is relatively prime to q, then 


K{qb,ab) = Y, ^-^P^E{q,ad 2 ). 

^ ^ cti 


did2=b 


9 






Proof. Since every divisor of qb can be written uniquely as a divisor of b times a divisor of g, 



which is equivalent to the statement of the lemma. □ 

The above quantities appear [6, Lemma 4] in a key part of Montgomery’s argument, which 
displays a bias in the values of Rq sampled on an arithmetic progression; 

Lemma 3.5 (Montgomery). There exists a positive real number c such that if 0 < a < q is 
a non-integral real number, then 

N 

Ro{nq + a) = C{q, a)N + O(^N exp(—cy^log N)'j 

n=l 

uniformly for q < ^. 

The proof of Lemma 3.5 is based on the following ancient identity due to Raabe [11] for 
a sum of B{x) over an arithmetic progression of points: for all real numbers x, 

J2B(x+fj=B{Jx). 

valid for all real x; the relevance to the problem at hand is due to the fact [6, Lemma 1] that 

d<x ^ ' 

To exploit Lemma 3.5, Montgomery chose a = g/4, so that for all divisors d of g the quan¬ 
tity B{a/d) equals ±|, with the sign depending only upon the residue class of d modulo 4. 
On the other hand, he chose g to be the product of many primes congruent to 3 modulo 4; 
for any divisor d of this g, the residue class of d (mod 4) depends only on the number of 
prime factors of d. With these choices, the sign of B{a/d) correlates exactly with the sign of 
fi{d), making the quantity K{q, a) large in absolute value. Choosing a = 3g/4 instead again 
makes K{q,a) large but with the opposite sign. 

One of our key observations is that we may work modulo a suitably chosen prime P rather 
than working modulo 4. Instead of choosing g to be the product of many primes congruent 
to 3 (mod 4), we instead choose g to be the product of many primes that are quadratic 
nonresidues modulo P. The sign of B{a/d) will not be perfectly correlated with p(d), but 
there will be enough of a systematic bias in the signs of B{a/d) (due to the imperfect 
distribution of quadratic residues and nonresidues modulo P) that we can still force K{q, a) 
to be large in absolute value. If we choose the prime P carefully, we can even force all of the 
different K{qbi, abi) to be large in absolute value and have the same sign. 

We introduce the following function, whose oscillations we will want to establish. For the 
rest of this section, we use boldface variables such as a to indicate the dependence of various 
quantities on /c-tuples (oi,..., a^) of variables. 

Definition 3.6. Given real numbers oi,..., and bi,... ,bk, dehne 

Ga,b(a;) = aiRoibix) H-h akRo{bkx). 
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The starting point of our modification of Montgomery’s method is the following easy conse¬ 
quence of Lemma 3.5. 

Lemma 3.7. Let Oi,..., and &i,. .., 6^ > 0 be real numbers. There exists a positive real 
number c such that: if 0 < a < q is a real number such that none of the abj is an integer, 
then 

N k 

Ga,b(^g + a) = N'^ajC{qbj,abj) + Oa,b(A^exp(-c-\/logiV)) 

n=l j=l 

uniformly for max{6i, ... ,bk] < q < ^. 

We now state our oscillation result for Ga,b; after which we show how Theorem 1.5 is 
implied by it. Thereafter our only remaining goal will be to establish this proposition; 


Proposition 3.8. Let Oi,..., and bi,... ,bk be fixed positive integers. There exists a 
constant Ka,b > 0, and sequences QN,h < ^ and 0 < afj < QN,h defined for positive 

integers N, for which QN,h tends to infinity with N and 

N 

y^,Ga,h{nQN,h + ajf) = Ka.bA^Vkyio^ Oa,b(iV) 


y^Ga,b(y^QjV,b + ajv) = -«:a,bA^yioglogiV Oa,b(iV)- 

77=1 


Proof of Theorem 1.5 assuming Proposition 3.8. By Lemma 3.2, 

k k 

'^rjE^isjx) = x'^rjSjRo{sjx) + Or,s(a:exp(-c\/logx)). 
i=i i=i 

Let Ds be the least common denominator of the rational numbers Si,..., s^, and set Oj = 
DgVjSj and bj = DgSj. Replacing x by D^x, we obtain 

k 

(3.1) y^rjEAjsjDsx) = xG^^b{x) + Or,s{x exp(-c^log x)). 

i=i 

Let 0 < £ < Ka,b- If there existed an xq such that Ga,h{x) < (Ka,b — £)'\/log \og{Dsx) for all 
X > Xo, then we would have 

N 

Ga,h{nQN + «iv) < (^a,b “ £)N 

71=1 

= (Ka,b - £)N I^Vlog logiV Os ^ 0£,a,b(l) 

by the bounds 0 < aff < QN,h', for large N (and hence for large QN,h) this would contradict 
Proposition 3.8. Therefore no such xq can exist, in which case equation (3.1) implies that 
there are arbitrarily large values of x for which 

k 

'^rjE^{sjDsx) > x{Ka,h - e:)yfog log^D^x) + Or,s {x exp(-cy/fog x)). 

1=1 



J\og\og{Ds{NQN + a^)) + 0{xo max Ga,b(^)) 

V l<i<xo 
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In other words, 


lim sup 


riE/^{six) + 


x^-oo Xyjlog log X 


^ ^a,b 

- lX'‘ 


and the analogous argument using the values G{nQ]\f + a^y) gives 

riE^isix)^ - \-rkEA{skx) ^ Ka.b 


lim inf 

X—)-oo 


x\/log logx 




completing the derivation of Theorem 1.5 from Proposition 3.8. 


□ 


Before addressing Proposition 3.8 directly, we record some preliminary facts about the 
distribution of primes in residue classes and the associated L-values. 

Lemma 3.9. Let P = 3 (mod 4) be a prime exceeding 3, and let xi{') = (p) denote the 
quadratic character modulo P. Then the class number h{—P) of the field Q(\/—P) equals 

H-P) = ^L{hXi) = -- 

TT Jr^ 

a=l 

Proof. These results are classical; see for example [?, chapter 6, equations (15) and (19)]. □ 


Lemma 3.10. Let P = 3 (mod 4) be a prime exceeding 3. If Xo denotes the principal 
character (mod P), then 


n 0+^) 

p<y ^ 


cp^/logy + Op{l), 


where 

(3.2) 


Cp 


P-1 

W h{-P)y/P 




1/2 


On the other hand, for any nonprincipal character x (mod P), 


n 0-f) 

p<y ^ ^ ^ 


<Cp 1. 


Proof. Again let Xi(-) = (p) denote the quadratic character (mod P). For any character 
X (mod P) we can write 


(3.3) n 


p<y 


Id 


x{p) 


P 


n 

p<y 


1 - 


x{p) 

P 


-1 


n 

p<y 


1 - 


x{p)xi{p) 


P 


n 

p<y 


1 - 


x^ip) 


pz 


The last product is absolutely convergent uniformly in P; indeed. 






<n 

n>y 


1 + 




<n 

n>y 


1 - 
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and so the last product equals n(^)=-i ~ x'^{p)p (l + ^(y)) particular is uni¬ 

formly bounded. When y is nonprincipal, we know [7, Theorem 4.11(d)] that 


n 

p<y 


x{p) 


P 


— L{1, x) + 


logy 


= L(l,x) 1 + 0 


\ogy 


since L{l,x) 7 ^ 0 [7, Theorem 4.9]. (Better error terms are available but are not relevant 
for us.) In particular, when neither x noi' XXi is principal, the first two products on the 
right-hand side of equation (3.3) are L(l,+) (l + Op and L(l, xxi)"^(l + Op(j^)), 

respectively; in particular, both are <Cp 1. This estimate establishes the lemma when y is 
neither principal nor equal to Xi- 

When X = Xi, the hrst and third factors are still bounded, while now the second factor 
actually diverges to 0, hence in particular is still bounded. Finally, when x = Xo is principal, 
the second factor converges to 1/L(l,xi) = \/P/vrh(—P) by Lemma 3.9, while the hrst 
factor on the left-hand side of equation (3.3) is 



Xo{p) \ 
P ) 


-1 


n 

p<y 

p^p 


1 - 


P 


-1 


P-1 

P 


(e^ logy) 




by Mertens’s formula [7, Theorem 2.7(e)]. This asymptotic evaluation of the right-hand side 
of equation (3.3) establishes the lemma when x is principal, indeed with the stronger error 
te™Op(y7I=)- ° 


We can now dehne the modulus appearing in the statement of Proposition 3.8, in 
terms of a companion prime Pb. 


Definition 3.11. Let Pb be the smallest prime satisfying Pb = —1 (mod 8bi...bk). Note 
that Pb = 7 (mod 8), and so —1 is a quadratic nonresidue (mod Pb) while 2 is a quadratic 
residue (mod Pb). Furthermore, if p is any odd prime dividing one of the bj, then by quadratic 
reciprocity [8, Theorem 3.1] we have, since Pb = 3 (mod 4) and Pb = —1 (mod p), 

Consequently, each prime dividing every bj is a quadratic residue (mod Pb). 

Now dehne, for any integer iV > 1, 


p<Cy/\og77 



where c is the constant from Lemma 3.7. The prime number theorem for arithmetic pro¬ 
gressions to a hxed modulus (see [7, Corollary 11.21]) tells us that logQAr^b ~ |c\/logiV; in 
particular, Qjv.b tends to inhnity with iV, and QN,h < when N is large enough. 

Note also that QN,h is squarefree and relatively prime to Pb and to each bj, since every 
prime p \ bj satishes (-^) = 1. Finally, note that any divisor d of QN,h has the convenient 
(and, for us, crucial) property that = p{d), since both quantities equal (—1 )#IpI'^I. 
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Lemma 3.12. For any 1 < b < P — 1, we have 

d=b (mod Pb) 

where Cp was defined in equation (3.2). 


Proof. From the orthogonality of the characters modulo Pbj 


V i 

, d 

d\QN,h 

d=b (mod Pb) 


' E E 


Pb — 1 ^ d 

X (mod Pb) d|Qjv,b 


E x(b) n (i 


Ph-i 


X (mod Pb) p|Qjv,b 


x{p) 


P 


Pb-1 


E x{b) n (i 

X (mod Pb) p<c^logN 


x{p) 


P 


Pb-1 


^—cp^^\og{c^/\ogN) + Op^{l) 


by Lemma 3.10. The statement of the lemma follows upon noting that log(c\/logiV) = 
iloglogiV + 0(l). □ 


Lemma 3.13. For any integer m that is not a multiple of Py^ and for any integer N > 3, 


K 



mQN,h \ 

) 




cp^K-Ph) 
V2(Pb - 1) 


V^log logiV + (9b(l). 


Remark. The exact value of the leading constant is not as important for us as the fact that 
its dependence on m is only in the term (^), so that the sign of the leading constant depends 
on whether m is a quadratic residue or nonresidue modulo Pb. 


Proof. By Dehnition 3.3 and the coincidence between the Legendre symbol and the Mobius 
function on divisors of QN,h (as noted at the end of Dehnition 3.11), we have 


-K 




1^/ mQN,h '\ 
d \ dPb / 



V i. 

d 

<i\QN,b 

mQjv,b/<i=“ (mod Pb) 
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The last congruence is equivalent to d being in the reduced residue class mQ^^h/o, (mod Pb), 
and so Lemma 3.12 applies: 


[ mq\ 


V 

mQN,h '\ 


Pv 


B 


a 




Pu-1 


-AglogivT — B — 

Bb Jv2(p^~iy ^ \p>, 

By the definition of the Bernoulli polynomial B, this sum equals 


Ob(l). 


Pb-i 


Pb 


Pb-i 


Ph 




P, 


Pb-i 


a=l ^ ^ ' a=l ^ ' a=l 

by Lemma 3.9, which completes the proof of the lemma. 


Pb 


□ 


We now have all the tools we need to establish Proposition 3.8 and hence our main theo¬ 
rems. 


Proof of Proposition 3.8. Given a sufficiently large integer iV, define two numbers = 
pP^Qn,^/P\>^ where the integers 1 < < Pb — 1 satisfy = iQv.b (mod Pb); note that 

none of the numbers a'^bj is an integer, since neither QN,h nor any of the bj is a multiple 
of the prime Pb. Since we confirmed in Definition 3.11 that Qv.b < we may invoke 

Lemma 3.7: 

N k 

y^Ga.b(nQjv.b + aff) = N^ajC{QN,hbj,affbj) + 0{Nexp{-c^/\og N)) 

n=l j=l 

= Nj2‘‘,K(QN.tbi,aM^ n (l“^) 

i=i plQiv.bL ^ ^ ^ 

As noted in Definition 3.11, each bj is relatively prime to QN,h- Thus by Lemma 3.4 with 
our choice of a^. 


N 


y^Ga,b(nQjv,b + <y%) = 0{N) 


n=l 


6 


TT^ 


PiQjV.b 

By Lemma 3.13, 

N 

^^Ga,h{pQN,h + 

n=l 


p 


-1 k 


p-L E“.np-k E 


i=l p\bj 


6 


P 


-1 


O^n) ~ 


X 


E 

did2=bj 


p\QN,h 

p{di) f frn^d' 


did2=bj 


— I k 


/i(di 

di 


K I QAr,b, 


m=^d2Qjv,b 

Pb 


o(A)+^A n (i-^j E“.n(i-^ 


t=l p\hj 


P 


-1 


di 


Pb 
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As noted in Definition 3.11, every prime dividing bj is a quadratic residue modulo Pb, which 
implies that (■^) = 1 always. Moreover, = Q]y,h (mod Pb), so the product of Legendre 

symbols (^) (^) equals 1; on the other hand, since —1 is a quadratic nonresidue modulo 
Pb, the product of Legendre symbols (^) (^^) equals —1. Consequently, 


N 


+ = ^^^ iVv/loglogiV J] (l-^) 

n=l ^2(Pb - 1) V P J 


— I k 


X 


p\bj 


^' di 

^ ' did2=bj 


The last sum equals X]d|fe ~ which when multiplied by the preceding product 

becomes Hpib (1 + addition, by the same argument as in the proof of Lemma 3.10, 


P\QN,h 


n (i-y 


-1 


hp =-l 


n L-L 


p 


-1 


VlogiV 


We therefore see that we have established the proposition with 


^a,b 


6 Cp^h{-Py,) 
vr2 V2(Pb - 1) 


n L 


=-i 




E«in(i+ 

i=l v\bj ^ 


P 


(3.4) 


e^/23^ h(-Pb)^/2 

7^5/2 _ 1 


hC =-l 




-1/2 k 


n l- 3 E%nL+; 


i=i p\hj 


P 


-1 


by the definition (3.2) of cp. 


□ 


Remark. The condition that all of the Oj are positive is used here only to ensure that this last 
line is non-zero. In the case where the aj may be arbitrary, and thus where the innermost 
sum in equation (3.4) may equal 0, it is possible that modifying QN,b and the above argument 
could prove the totient independence for any coefficients; however, we succeeded in getting 
this to work only when k < 3. 
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